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We study experimental signatures of tiie Ising nematic quantum piiase transition in d-wave su- 
perconductors, associated witli tlie cliange of lattice symmetry from tetragonal to orthorhombic in 
the superconducting state. The characteristic feature of this transition is that the ratio between 
the Fermi velocity vp and gap velocity flows to a maximally anisotropic fixed point, i.e. the 
renormalization group fixed point is situated at (va/vf)* ~ 0. Our main point is that the logarith- 
mic approach to this fixed point has visible signatures in the thermal transport. The analysis of 
the quasiparticle contribution to the thermal transport is carried out in the framework of a kinetic 
approach, which shows that the thermal conductivity is enhanced near the nematic critical point. 
Another aspect of our study is the interplay of dilute disorder and electronic interactions in the 
measured thermal transport coefficients. 

PACS numbers: 



I. INTRODUCTION 

Electronic nematic phases were discussed in the con- 
text of doped Mott insulators^, and have by now been 
experimentally observed in a number of systems. They 
have been observed in semiconductor heterostructures^, 
in the bulk transition metal oxide Sr3Ru207^, as well as 
in YBCO. In all these cases the thermal transition to the 
nematic phase seems to be second order, whereas the na- 
ture of the quantum phase transition is much less clear. 
In the case of Sr3Ru207 it seems to be first order^.. 

Now there is evidence that a nodal nematic phase oc- 
curs in at least some of the underdoped cuprate supercon- 
ductors. In the nematic phase, the square lattice sym- 
metry is broken down to tetragonal symmetry, a con- 
sequence of the instability of the interacting electronic 
system to partial stripe-like order. The best evidence 
for this comes from measurements of strongly temper- 
ature dependent transport anisotropics in underdoped 
YBa2Cu306-i-5^ and from neutron scattering experiments 
in underdoped YBa2Cu306.43^. 

On the theoretical side, Vojta et al.& analyzed possible 
quantum phase transitions in d-wave superconductors in 
the framework of a renormalization group (RG) analy- 
sis. An initial RG analysis found runaway flow for the 
nematic ordering instability at zero temperature based 
on an expansion in 3 — d, where d denotes the spatial di- 
mensionality. In a recent work, Kim et al^ found in the 
framework of a large-N analysis the existence of a second 
order transition to a nematic phase. Based on an RG 
analysis in the lavge-Nf framework, where Nf denotes 
the number of electronic spin components (the physical 
case corresponding to Nf = 2), Huh* et al. confirmed 
the existence of such a second order transition. The au- 
thors found an RG fixed point at order 1/iV/ describing 
a second-order quantum phase transition associated with 
the onset of long-range nematic order. The scaling prop- 
erties near the fixed point are very peculiar. There is a 
dangerously irrelevant parameter v/^/vp, where ua and 



vp are the velocities of the nodal fermions parallel and 
perpendicular to the Fermi surface, which controls the 
fixed point. The fixed point lies at "infinite anisotropy", 
i.e. (va/vf)* — 0, which has to be contrasted from 
the other fixed points found for other competing orders, 
which are relativistically invarianl*^. In order to calculate 
physical quantities one has to use a fully two-dimensional 
theory, since the flow of the anisotropy to zero is logarith- 
mically slow as a function of the relevant energy scale, as 
is described by Huh et aZ..~ 

In the present paper we study thermal transport prop- 
erties at the nematic to isotropic quantum phase tran- 
sition (QPT) deep within the d-wave superconducting 
phase of a quasi two dimensional tetragonal crystal^ in 
the framework of the Boltzmann equation. The main re- 
sult of our analysis is the logarithmic enhancement of the 
thermal conductivity upon lowering the temperature, see 
Eq. p.33p and its disucssion in the result section together 
with Fig. H 

The paper is organized as follows. We start with a 
general review of the model for the second order ne- 
matic phase transition and its properties in a few limiting 
cases in Sec. [Til This closely follows the presentation in 
Refsjii^. A review of the properties of the model under 
renormalization group transformations is given in Ap- 
pendixIXl We proceed with the deflnition of the heat cur- 
rent operator of the electronic quasiparticles in Sec. lIII Al 
as we use it in the present work. In Sec. IIII Bl we in- 
troduce the Boltzmann equation framework and explain 
the variational ansatz used to solve it. In Sec. IIII CI we 
review a Boltzmann equation analysis of a disordered 
d-wave superconductor without a nematic mode, which 
makes connection to existing results^SiSS,. Then we an- 
alyze the full problem taking into account elastic scat- 
tering from dilute impurities and inelastic scattering due 
to the nematic mode in Sec. IIII Dl The numerical solu- 
tion of the Boltzmann equation is presented in the results 
section. Sec. IIV Bl and several different situations are dis- 
cussed, such as the interplay of inelastic and elastic scat- 
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tering. Finally, in Sec.|V]we conclude and comment on 
possible experimental implications of our analysis. Ap- 
pendix [C] provides some additional information on trans- 
port in clean systems and complements the discussion in 
Sec. nil Al Appendix |B] further analyzes the problem in 
the framework of a Kubo formalism, and points out a 
couple of difficulties. This only serves to complement the 
picture and may serve as a starting point for more ambi- 
tious approaches to this complicated transport problem. 



II. THE MODEL 



where e^^-a is the antisymmetric tensor. 

The electronic Hamiltonian in Eq. (|2.ip has a Dirac 
Hamiltonian structure and is readily diagonalized with 
a standard Bogoliubov transformation. This is achieved 
by a transformation according to 



7Lk = "k/Lk-fk/; 



3-cr-k ' 



73-a-k ^ "k/3_^_k + ^^k/i^k 



(2.5) 



and respectively for nodes 2 and 4; itk and Wk are the 
coherence factors^. The corresponding unitary matrix is 
given by 



The model under consideration throughout this work 
has been discussed in the literature in great detaitfii^. 
Therefore we will only repeat the key features. The rel- 
evant low-energy description of the electronic system in 
a two-dimensional d-wave superconductor with a pure 
dg.2_y2 pairing symmetry is given by the following BCS- 
type Hamiltonian 
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where the dispersion has been linearized around the four 



nodal points, 
velocities v^. 



The Fermi velocities v 
are defined as 



vi and the gap 
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where the energies for the Bogoliubov quasiparticles are 
given by 
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In our later discussion of the thermal transport properties 
in terms of the Boltzmann equation it is mandatory to 
work in the basis of these quasiparticles^^'^^. 

In the theory considered^i^ii^ an Ising-symmetric ne- 
matic order parameter couples to the nodal fcrmions. 
The corresponding interaction term has the form of an 
additional s-wave order parameter, whose condensation 
has the effect of breaking the four-fold rotation symme- 
try in k-space down to a two-fold one^. Consequently, 
it assumes the following form 



and 



with V = — ^ 



(2.3) 



where the vectors denote the location of the nodal 
points in the Brillouin zone in clock-wise direction start- 
ing with Ki lying at (7r/2,7r/2). Furthermore, for later 
convenience, we introduced the anisotropy parameter 
iJ ~ vf/v/\ in the above equations. This parameter plays 
a vital role since it is a direct measure of the velocity 
anisotropy and has a nontrivial flow under the renormal- 
ization group transformation^. The main results of this 
RG analysis are subsummed in Appendix |^ Further- 
more, it is notationally very convenient to introduce the 
nodal fermions, called fia, /2(t, /so-, and fia living at 
the respective nodes in k-space. The index a denotes a 
spin degree of freedom, i.e. a =1,1. This can be ex- 
tended to allow for a large-N treatment by generalizing 
a = l,...,Nf. Our notation in the following will how- 
ever stick to the SU(2) version. The Nambu spinors are 
composed of the nodal fermions in the following form 



Nf 

E 



(fxdrcj) 



la 



2a 



,(2.8) 



where </> denotes the Ising-type order parameter, whose 
action can be obtained by integrating out the fermions. 
Kim et ali^ showed that the leading quantum fluctuations 
at the large-N level after a quadratic expansion around 
the saddle point lead to a non-analytic form of the effec- 
tive bosonic theory. This calculation is straightforward 
but tediousi^iii, relying on an application of Feynman 
parameters. Retaining only the terms which are relevant 
at low energies in the RG sense (assuming the mass of 
the bosonic action is tuned to criticality) , the effective 
bosonic action assumes the following form 
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where 
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(2.10) 



3 



This implies that the effective propagator of the bosonic 
modes is given by 



III. THE KINETIC APPROACH TO HEAT 
CONDUCTIVITY 




4ip) - Kp) 



4ip) 



(2.11) 



Higher order terms are as usual in large-N theories down 
in powers of 1/Nf. At this point it is instructive to dis- 
cuss the above propagator in two extreme cases, namely 
for TJ = 1, i.e. isotropic velocities, and for w ^ 1. 

For V — 1 the corresponding bosonic propagator is sim- 
ply given by 



Dt 



■e2(p) 



(2.12) 



where the energies ei and £2 are now identical due to 
the perfect isotropy. It is important to note that the 
present propagator implies the existence of quasiparticle- 
like peaks sitting on top of a continuum in the spectral 
function^. The infrared RG fixed point of the theory lies 
at infinite anisotropy. In the limit of extreme anisotropy, 
i.e. V ^ 1, the propagator reads 



^(V^^ + (vi.-p)' + V'^?. + (v|-p)') (2-13) 



It was pointed out before that the presence of the ne- 
matic order parameter and its interaction with the Bo- 
goliubov quasiparticles leads to a finite lifetime for the 
nodal quasiparticles^.. 

In order to discuss transport in a realistic setting we 
will also add some dilute disorder. On a Hamiltonian 
level this reads 



Hdis — 



(2^(2^ 



where Vj^^, stands for the scattering matrix element of 
an electron living around node i with momentum k into 
an electron living around node j with momentum k'. In 
contrast to Refs. we assume a source of isotropic 



scattering, which implies that V, 



kk' 



V^j., — u. However, 



all the results of the following calculations can in prin- 
ciple be extended to account for the more generic case 
including different scattering strenghts. As a remark it 
is important to note that our results will not recover the 
unitary limit. 



A. The quasiparticle expression for the heat 
current 



In our present paper we will concentrate on the con- 
tribution carried solely by the nodal quasiparticles. The 
heat current carried by the Bogoliubov quasiparticles is 
readily given by the intuitive expression 



N, 



d^k „ , de 
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(3.1) 



This expression will constitute the starting point for our 
Boltzmann transport equation analysis in Sec. IIII Bl 

For a general superconductor, irrespective of the pair- 
ing symmetry, we can refer the reader to Refs J^i^SiiSii^SiiH 
for a proper derivation of the appropriate operator. 

For an outline of a generic expression for the heat cur- 
rent in an interacting system we have added some re- 
marks in Appendix iBl which also addresses the question 
of the role of the bosonic fluctuations. 



B. The Boltzmann equation 

Within this section we access the thermal transport 
properties using the semiclassical Boltzmann-equation 
approach. This approach has proven to be a powerful 
tool to compute transport properties of quantum critical 



systems 
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In the following we assume that 



a quasiparticle description applies^. The controlling pa- 
rameter of the Born approximation is given by 1/Nf. 

The central object in Boltzmann transport theory is 
the distribution matrix of the quasiparticles. In our case 
those are the Bogoliubov quasiparticles. We introduce a 
distribution function of quasiparticles of the form 



/:(k,t) = (7t(k,t)7,,(k,i)) 



(3.2) 



For all our following considerations it is important to 
assume that the Bogoliubov particles constitute reason- 
ably sharp quasiparticles. In equilibrium, i.e. in the ab- 
sence of external perturbations (such as an applied volt- 
age, temperature gradient, ...), the distribution function 
is given by familiar Fermi-Dirac distribution 



/;(k,t) = n0(6,(k)) 



1 



1 



(3.3) 



where ei(k) is given in Eq. (|2.7p . 

In order to deal with all the Bogoliubov quasiparticles 
on equal footing, we introduce a local basis, where the 
node is again parametrized by i. We assume that we 
apply a temperature gradient across the system, such 
that the temperature at position r is given by T(r) = T-\- 
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r • VT. The Boltzmann equation assumes the following 
schematic form 



(3.4) 



and since we are only interested in the time-independent 
solution we arrive at the following simple equation 



X,(k)=J^,ii(k) 



(3.5) 



We will now carefully disentangle the different terms in 
this expression. The driving term Xi assumes the generic 
form (we drop the spin-index for reasons of simplicity; it 
will trivially be accounted for by a factor of 2 in the end) 



X(k) ^^'W^^, 
dedk) _ (vl.(v5..k)+vX(vX-k)) 



dk 



(k)n'^(e,(k))(l-n"/e.(k))), 

(3.6) 



e.(k) 



The form of Eq. p.6p is markedly different from an 
isotropic system. In an isotropic system it is always pos- 



sible to choose a basis such that one can formulate the 
problem in terms of the angle enclosed between VT and 
k, i.e., X, oc 4|VT|fccos(Z(Vr,k))^>25. This allows 
to find a much simpler solution than in our case, where 
the angular dependence has to be taken seriously. Our 
strategy to properly account for the anisotropy will be 
to consider two independent equations for the particles 
moving parallel and perpendicular to the nodal points. 

Furthermore ZcoU denotes the so-called collision term 
or integral. The collision integral Z^^^jj for a particle liv- 
ing at node i will in general be composed of two sources 
of relaxation, namely the contribution due to disorder, 
henceforth called Z^jg, and the contribution due to the 
inelastic scattering, i.e., 



22„,i(k)=2SJk)+2^,i(k) 



(3.7) 



The following sections are devoted to the understanding 
of the single and combined effect of those two scattering 
mechanisms. The specific form of the two contributions 
is easily obtained and reads 



X^i, = 27rn;n.pu'Yl J ^ {u^i.r.U-^') {U,,.'r.Ur^') ^^k) - e,{k')) [fik) - /^(k')] 



S {lj + uj' + €j{k'))lmD{uj' ,k') x 



X [r(k)(l + Pik + k')) + n(k') (r(k) - Pik + k'))] 

I 



(3.8) 



The set of equations Eq. jSZ]), dSJ]), is an- 

alyzed in the remainder of the paper. At this point we 
will perform a set of manipulations on the above expres- 
sions to bring them into a more tractable form. It proves 
convenient to introduce the following simplifications at 
this point. Following Eq. (|3.10p we can use the local 
basis at any nodal point, which allows to rescale mo- 
menta to yield an isotropic coordinate system, where the 
anisotropy is absorbed in the measure of the integral. 
We furthermore rescale momenta by a factor of T which 
consequently eliminates all factors oi vp, va, and T ex- 
cept for global prefactors if scattering processes within 
one node or across the Brillouin zone are considered. If 
scattering processes mix adjacent nodes, however, this is 
not fully possible and factors of v prevail. The rescal- 
ing also allows the use of polar coordinates and implies 
that energies are now given by = ±kT. This allows to 
rewrite the scattering term locally as 

VT 

X, = -(v^e.+vkA.)fc— n?(fc)n5.(-fc) 

= Xf + X^ (3.9) 



where 

m = =cos(0 + (z-i)|) , 

A,{0) = sin(0+(z-l)|) , (3.10) 

rendering all expressions dimensionless. The driving 
term splits into two parts, one characterizing parallel 
and one perpendicular electrons which is consistent with 
treatments using the Kubo formula^^. 

The structure of the driving term motivates an ansatz 
for the solution of the Boltzmann transport equation 
given by 

V* VT 
V* VT 

+ vXA.(0)^^4(fc)nO(-fc)vI/X(fc,0,«), 

(3.11) 

where the superscript i of the functions 5'^(fc,0,U) and 
\I'^(fc, d,v), respectively, accounts for the fact that we can 
define a local basis for every Dirac point itself. Again, it is 
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worthwhile to contrast this expression from an isotropic 
system. In such a system we would not need to dis- 
tinguish the two different nodal directions and we could 
simply choose an ansatz of the form ^fdlcDyT • k due to 
the spherical symmetry of the problem. 

A few words on the symmetry of the obtained expres- 
sions are in order here, since this will help to simplify 
life a lot in the following. The symmetry of the driving 
term under exchange of vis-a-vis nodes, i.e. Xi{k,9) = 
Xi+2{k,6), enforces Sfi{k,9) — 6fi+2{k,9). Furthermore 
we can deduce from the symmetries of the driving term 
that 

(3.12) 

Following Eq. p.ip we can derive an expression for the 
energy current carried by the Bogoliubov particles of the 
form (note that we have Nf = 2 in this expression) 

iE ^ ^izf ^fcMvK.+vkA.]<5/.(fc,0) 



^ f dn^dk 



vpVA ^ J (27r) 

z— 1 



(3.13) 

Taking into account all the aforementioned symmetries 
of the problem we finally arrive at a relatively simple 
expression for the thermal current carried by the quasi- 
particles under an applied thermal gradient across the 
sample 



+ Va 



^krde)n"f{k)n",{-k)^Uk,0,v) 
dCl]s_dk 



k^^l{e)n){k)n]{~k)^\{k,9,v) , 

(3.14) 

where the terms mixing vp and wa vanish according to 
Eq. dSUl), due to 



de/^,{e)£,,{e)^'p/^{k,e,v) = o . (3.15) 



This is analogous to the vanishing of mixed terms in the 
treatment by Durst and Lee^^. We furthermore intro- 
duced the following abbreviation: 

4T2 



V 



_F/A 



VpVA 

(3.16) 

Since we are interested in the transport coefficient k^x 
we can give its generic expression as 

Kx. = -4T^v I ^^k^eiie)n%k)n'f.i-k)^Uk,0,v) 

- ^/ ^k'Al{e)n'}{k)n'}i~k)^l{k,e,v) . 

(3.17) 



C. Thermal transport in a disordered d-wave 
superconductor 

Within this subsection we investigate the thermal 
transport in a system where the only relaxation mech- 
anism is provided by dilute disorder. This has also been 
discussed in the thesis by Paaske^^ and simply serves as 
a reference. We introduce a parameter measuring the 
strength of impurity scattering 



2n 



where 



a 



imp 



(3.18) 



(3.19) 



Expanding the collision term due to disorder to linear 
order in the deviation from equilibrium leaves us with 
the following expression 



ikT f dQh 



2n 



4 

■ ^T+(0,^) SMk,e) 



<5/,(fc,(/.) 



(3.20) 



where we introduced the short-hand notation 



1 + <.(^)^,(0) - kA,(0)A,(0) , (3.21) 



with ^4(6*) and Ai{9) defined in Eq. ((XTU)) and k = ±. 
This factor simply accounts for the usual coherence fac- 
tors2.. We are thus left with the task of solving the equa- 
tion 



x,ik,0)^T^,,{k,e) . 

Plugging in the ansatz (|3.1ip we see that 



X^i,(fc,0) = akTSUk.O) 



akTSf,{k,e) 



ak 



E 



(3.22) 



-Sfjik,^) 



(3.23) 



where the second term in the first line vanishes due to 
the isotropic nature of the impurity scattering, which we 
assumed for reasons of simplicity. In the more generic 
case--- this is no longer true. In our problem, however, 
the second term will always vanish also in the presence 
of the inelastic scattering by virtue of Eq. (|3.12p . 
We can now find a simple solution given by 



1 



^Uk,O,v)^^\{k,0,v) = - 



a 



(3.24) 



Using Eq. (I3.17|) it is straightforward to arrive at an ex- 
pression for the thermal transport coefficient, given by 



"""" " — P ^''^ + ""^^ 

^ ^ ' ''imp ^ 



(3.25) 
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which is the central result of this section. Paaske^S. 
showed that this discussion can easily be extended to in- 
clude the case of anisotropic scattering, which means we 
are formally allowing for three different scattering matrix 
elements: intranodal scattering, which we denote V^^,, 
scattering between adjacent nodes called V^^,, and V^^, 
for scattering across the Brillouin zone. In this more gen- 
eral case the thermal conductivity reads 



ttT 



12n. 



imp "'O 



i + s 



(3.26) 



where 



2 {V^i^. 



2ul 



(3.27) 



This expression nicely reduces to Eq. ()3.25p for isotropic 
scattering, i.e., V^^ = V^l - V^^,. 

In the following sections we will concentrate on 
isotropic scattering to study the interplay of disorder 
scattering and scattering from the nematic order param- 
eter for reasons of simplified analysis. Another comment 
on this result is in order here. It has been shown that in 
the context of unconventional superconductor a faithful 
description of disorder requires to consider unitary scat- 
terers, see Ref.-iS. This means that a finite concentra- 
tion of non-magnetic impurities induces a finite density 
of states. The so-called universal conductivity obtains in 
the limit when temperature T is much smaller than the 
impurity bandwidth 7, i.e. T <^ 7, see Ref. [2^ The 
universal heat conduct ivity^i^ has been shown to be in- 
dependent of details of the disorder distribution, which 
means that the impurity concentration drops out in the 
final expression for the heat conductivity. Experimental 
evidence also points towards the existence of a univer- 
sal conductivity in the limit T — > independent of the 
scattering strengt h^^i^^i^^ . 

Here, the result depends upon the impurity concen- 
tration explicitly. Our calculation, however, does not 
address the unitary limit but the Born limit, such that 
this discrepancy does not constitute a problem. In or- 
der to make contact with the unitary scattering limit. 



the above treatment must be (in the spirit of a quantum 
Boltzmann equation^) supplemented by a field renor- 
malization stemming from the real part of the self-energy 
(this is neglected in the above Boltzmann equation, but 
can be incorporated in a straightforward manner^'''^^), 
which cancels the explicit dependence upon the impurity 
scattering. This establishes the equivalence of the Kubo- 
formula calculations employing a self consistent Born ap- 
proximatioi*i^, where the impurity scattering induces a 
finite density of states at the Fermi level. 

D. Thermal transport at the nematic transition 

Within this section we consider the electronic scatter- 
ing off the nematic mode. The central approximation in 
this section is to assume the bosonic sector to be in equi- 
librium, which implies that the bosonic system is not 
excited by the temperature gradient and relaxes on a 
time scale faster than the time-scale associated with the 
electronic quasiparticles. One argument in favor of this 
point of view is to consider a situation in which the tem- 
perature gradient is applied along the vp direction of the 
electrons at the nodal point denoted 1. Under the formal 
assumption that the rest of the electronic quasiparticles 
is accelerated by the driving field in their ua direction, 
they are in the limit of large anisotropy subdominant in 
their contribution to the thermal transport. The propa- 
gator of the bosonic mode is linked to contributions from 
all the fermionic nodes, and so its response to the ther- 
mal gradient is suppressed by a factor of 1 /N, where N is 
the number of nodes. An analogous discussion has been 
carried out in the study of transport in bosonic theories 
in the large-N limitW,. 

Another argument in favor of this approximation 
comes from rewriting the problem integrating out the 
bosons, which leads to a purely fermionic problem, see 
Appendix [C] In the fermionic language one can clearly 
see that a scattering involving fermions at adjacent nodes 
can relax the thermal current, whereas scattering be- 
tween electrons at the same node or nodes across the Bril- 
louin zone cannot. This argument, however, only works 
for large anisotropy, i.e. zJ 3> 1. 

We first elaborate on the collision integral stemming 
from the scattering of the nodal fermions from the ne- 
matic order parameter fluctuations. It assumes the fol- 
lowing generic form 



= I ^dk'k%-,{e, c^)D';^t-J' {{uBik - k') + n){-k')) Sf\{k, 9) - {usik k') + n){k)) SMk', 0)) 

+ / ^dk'k%-,(e,cb)D';^ti'^' [{nB(k + k')+n'}{k'^^^ , 

(3.28) 



which is a generalization of the expressions shown in Refs. accounting for the coherence factors and the 
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different nodes. Furthermore, we exploit the fact that 
6fi{k,9) = Sfi+2{k,d) and introduced D^^'J^ being 7 
times the imaginary part of the retarded Green's func- 
tion of the bosonic modes, which was introduced in its 
imaginary frequency form in Eq. (j2.1ip . The full problem 
is analytically not tractable and has to be solved numeri- 
cally. We thus use a variational approach^i^i^i^, which 
allows to determine a bound for the conductivity. The 
interaction parameter A drops out of the problem exactly 
which simply reflects the fact that the present perturba- 



tion theory is not controlled in the smallness of A, but 
in the smallness of l/Nf. It is interesting to note that 
the functions ^F/A{k, 0,v) acquire a true angular depen- 
dence in contrast to the pure isotropic impurity scatter- 
ing problem. The presence of the ncmatic mode leads to 
a non-trivial renormalization of the velocity-parameter T7, 
see Ref. [Sj, which will be taken into account later. 

We set up the variational problem defining the appro- 
priate matrix elements 



Q = Qf + Qa=X~M~V 

X = -— ^ j ^cose^^p{k,e,v)x!{k,6,v) + -—^ j ^M\{k,e,v)x^ 

AT r kdkdf) dth ~ ,, , 9 

M = — j ^^^dk'k'T~^{9,(p)D^t-^'; n]{k)n'>f{^^^ 

^ %j '^^^dk'kXr:(e,cf>)D;^t-J^'n°f(k)n'}^^^ 

AT /■ kdkdf) d(h 9 
+ j^J j^^dk'kXr:^,{e,4,)D^ti'^n'}{k)n}{-k){nB{k + k'^ (cos 0, IJ) - cos 

aT /•fcg^^^0(fc)„0(_fc)(^^^^^.^(fc^^^-))2^^ j h^kn){k)n]{-k)(sine^\{k,e,v)Y , (3.29) 



^ 2 



where X encodes the driving term, TV denotes the scatter- 
ing from the nematic mode, and 2? is simply the scatter- 
ing from impurities. The above form is relatively simple 
and its simplicity stems from the fact that there are a 
couple of different symmetries which can be exploited. 
In the following calculation we will assume Nf = 2. In 
this framework it is possible to obtain the Boltzmann 
transport equation by demanding a maximization of the 
functional in the following sense 



dcosO-^'pik.e.v) dsme^\{k,9,v) ' ^ • ; 

By virtue of Eq. p.lSp the contribution associated with 
the direction perpendicular to the Fermi surface and the 
one associated with the parallel (gap) direction decou- 
ple nicely, and we are left with two equations that wc 
can solve independently. Furthermore, comparing the 
expression for X in Eq. (I3.29P and the expression for 
the transport coefficient Kxx given in Eq. (|3.17p we see 
that they are identical up to prefactors. This allows to 
determine k^x knowing X. In order to solve the above 
integral equation we can now make a variational ansatz 
and maximize the above functional with respect to a set 
of expansion coefficients, allowing to determine a lower 
bound on the contribution'*^. It seems sensible to choose 



the following ansatz 



^'p/^{k,9,v) = -~-aF/A{v)^bF/A{v)cos2e , (3.31) 



which again amounts to dealing with fermions perpen- 
dicular and parallel to the Fermi surface, separately. In 
the above expression Qf/a bF/A £^re the variational 
parameters, which have to be determined as to minimize 
the functional. 

Of course, the whole series cos 2nd with n being an in- 
teger is allowed for symmetry reasons. However, it is easy 
to check that only n — 1 contributes to the heat current, 
which is why wc concentrate on this mode. We choose 
the mode with no k-dependence, since this is the mode 
associated with energy conservation in a clean system, 
see Appendix [C] This mode in a clean system with a 
fully relativistic Hamiltonian is not relaxed and leads to 
an infinite thermal conductivity. One can of course im- 
prove upon the approximation by including more modes, 
which is however beyond the scope of this work, since we 
are mainly interested in qualitative features. It is impor- 
tant to note that the expressions for X and V can be 
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calculated analytically, yielding 

TT 



a&{viT)) + 



X 
V 



24 
aTn 

Is' 



dp + OA 



2 



apbp + flA&A 

(3.32) 



The thermal conductivity will finally be given by the sim- 
ple expression 



6 

TT 1 

6 V 



bpiv) 



apiv) 

, bA{v) 



(3.33) 



At this stage it is worthwhile stressing the fact that the 
anisotropy parameter scales at criticality like 



v(T) = log 



(3.34) 



It is an interesting question to analyze the role of the 
"vertex" -corrections, which are encoded in the coeffi- 
cients OA, ap, 6a 7 and bp. Having a closer look at 
Eq. p.33p reveals that for large initial anisotropy the 
second term will be subdominant. Asymptotically the 
heat conductivity is thus given by 



ap{v) 



bp{v) 



(3.35) 



A closer inspection of the numerics reveals that the 




Figure 1; (Color online) The uppermost curve (solid/blue) 



shows the temperature behavior of V{T) = ^Itq) • 
lower curve (short dashes/red) shows ratio A{T) 



The 



plotted with an initial value of vf/va ~ 



20 at a temperature Tq. 

fiow of the velocity anisotropy and the flow of the "ver- 
tex" (the expression in the brackets) have the same log- 
arithmic dependence upon temperature. This is shown 
in Fig. [U where V{T) = is compared to A{T) — 



This implies an asymptotic behavior 



of the heat conductivity according to 



oc log — 



(3.36) 



The role of the vertex is also highlighted in the subse- 
quent discussion and in Fig. [2] 



IV. THERMAL CONDUCTIVITY AT THE 
NEMATIC PHASE TRANSITION 

We now turn our attention to the numerical solution 
of Eq. (I3.33|) in various situations. 

In a first part, Sec. IIV A[ we will analyze a clean sys- 
tem, in which the only current relaxation stems from in- 
elastic scattering of the electronic quasiparticles from the 
effective bosonic mode. In a second step. Sec. IIVBI we 
will add dilute disorder to the problem and consider the 
full problem. 

It was shown that the nematic phase transition is de- 
scribed by a fixed point located at infinite anisotropy, i. 
e. vp/vA — > oo^ . This fixed point is approached in a 
logarithmic manner upon lowering the temperature. This 
implies that in order to do a realistic calculation of the 
thermal conductivity one has to take into account the log- 
arithmic flow of the velocity-anisotropy, which will lead 
to a logarithmic enhancement of the thermal conductiv- 
ity upon lowering the temperature. 

The RG flow equations of the velocities which need to 
be integrated to achieve this are reviewed in Appendix [XI 

All the plots presented in this section originate from 
a combination of a numerical solution of Eq. p.33p and 
the flow equations presented in Appendix [Xl 



A. Quasiparticle thermal conductivity in a clean 
system 

In this section we consider a clean system, in which 
the only source of thermal current relaxation stems from 
inelastic scattering. Experimental evidence suggests a 
value of vp/vA ~ 20. Our reference point in the fol- 
lowing plot is thus given by this ratio; since we are only 
interested in the qualitative features the initial tempera- 
ture of the integrated flow is chosen arbitrarily and called 
Tq. From Eq. (|3.33p it is obvious, that there are differ- 
ent contributions, which we aim to disentangle. One is 
the obvious dependence on v and thus on its asymptot- 
ical flow (see Eq. (|3.34p ). Another one comes from the 
implicit dependence of the parameters of, a and ap^A on 
V and thus on T. In order to disentangle this, we have 
Fig. [2] contain a comparison of three curves in the clean 
limit. The uppermost curve is a full numerical solution 
of Eq. (|3.33p , where all parameters and their respective 
RG flow are fully incorporated, i.e. the full dependence 
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of a,b on v, and thus implicitly on T is taken into ac- 
count. The middle curve serves to contrast this from a 
situation, in which a, b would only depend upon the ini- 
tial value of V, but their flow is not taken into account. 
This curve thus simply reflects the flow of the prefactors 
in Eq. (|3.33p alone and thus asymptotically mimics the 
flow of the velocity anisotropy alone, see Eq. p.34p . The 
third curve simple serves as a reference point showing 
an unrenormalized flat curve, denoting the value of the 
thermal conductivity obtained at Tq with vf/va = 20. 



140 
135 
130 
125 



a=0 
a=0.5 
tt=1.0 
a=2.0 



0.05 0.10 0.15 0.20 0.25 0.30 0.35 T„ 



flow of anisotropy 



unrenormalized curve 



.10 -\ 




0.05 0.10 0.15 0.20 0.25 0.30 0.35 r„ 

Figure 2: (Color online) The uppermost curve (solid/blue) 
shows the full numerical solution of Eq. (|3.33|l . taking 
into account the temperature renormalization effect of the 
anisotropy ratio. The middle curve (short dashes/red) 
shows the universal limit conductivity with running couplings, 
whereas the lowest curve (long dashes/green) simply serves as 
a reference curve. All plots employ an initial anisotropy ratio 
vf/va ~ 20 at a temperature Tq. 

In the next section we will additionally consider the in- 
terplay between disorder and interaction with the bosonic 
mode. 



B. Quasiparticle thermal conductivity in a 
disordered system 

As we showed in the discussion of the thermal conduc- 
tivity in a disordered d-wave superconductor, the impu- 
rity scattering strength does not vanish if the field renor- 
malization is not taken into account, see Sec. IIII Cl This 
of course has consequences if one considers the problem 
taking into account both elastic and inelastic scattering, 
as shown in Fig. [31 We see that for different disorder 
levels, parametrized by a, the curve has different offsets, 
being maximal in the clean system, as one would expect. 
On top of this, by the definition of a given in Eq. (|3.f 8p . 
it is obvious that the disorder strength is a flowing pa- 
rameter under the renormalization group transformation, 
since it explicitely depends on the velocities. This effect 
is further illustrated in Fig. 21 where full solutions of 
Eq. (|3.33p are plotted. The difference is that the upper 



Figure 3: (Color online) This plot shows different impurity 
strengths in increasing order from top to bottom with the 
uppermost curve being for the clean system {a = 0, a = 
0.05, a — 0.1, Q = 0.2), i.e. Kxx/T, for an initial anisotropy 
ratio vf/va = 20 at an arbitrary temperature To for different 
disorder levels. 



curve takes into account the disorder renormalization, 
whereas the lower curve, which is plotted for reference, 
neglects this effect. 



renonnalized disorder 



\ 



Figure 4: (Color online) This plot shows the effect of disorder 
renormalization (initial value a — 0.01) due to the RG flow 
of the velocities. The upper curve (solid/blue) shows the full 
solution of Eq. (|3.33|l taking into account the disorder renor- 
malization, whereas the lower curve (dashed/red) is plotted at 
fixed disorder, but everything else renormalized. Again, i.e. 
i^xx/T is plotted for an initial anisotropy ratio vf/va ~ 20 
at an arbitrary temperature Tq. 



V. CONCLUSION 

In this paper we have shown that at the nematic tran- 
sition there is an enhancement of the thermal conduc- 
tivity due to renormalization effects towards the fully 
anisotropic RG fixed point at (va/vf)* ~ 0. The central 
result of this paper is thus given by Eq. (|3.33p and shown 
in Fig. [2l We addressed this experimentally relevant is- 
sue in the framework of the Boltzmann equation for the 
fermions, arguing that thermal transport of the collective 
bosonic excitation is suppressed in the large Nf limit^. 
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Even though our approach strictly speaking is only jus- 
tified under the conditions explained in Sec. IIII Di we 
believe that our treatment correctly captures the loga- 
rithmic enhancement as the central quantitative signa- 
ture of the nematic transition. The full treatment in- 
cluding the bosonic mode is beyond the scope of this 
paper and requires a full discussion in the framework of 
coupled Keldysh equations or a treatment in the frame- 
work of the Kubo formula or equivalently the memory 
matrices^. 
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Appendix A: RG EQUATIONS 

The full details of the derivation of the renormaliza- 
tion group equations were presented elsewhere^, so we 
only repeat the crucial formulae. The methodology of 
the applied RG approach differs slightly from the more 
standard hard cutoff scheme, since a soft cutoff is more 
favorable in order to deal with the anisotropic veloci- 
ties. Under the renormalization group transformation 
the fermionic velocities vf and va modify according to 



dvp 



= (Cl - C2)VF 



and 



(Al) 



(A2) 



This implies that the ratio of the two velocities scales like 

d{vA/vF) 



M 



{C2-C3){va/vf) 



(A3) 



where the functions Ci, C2, and C3 are given by 



Ci{v) 



2v- 



— 1 /•oo /•27r 



V J-00 Jo 



dx / def+--{x,0,v)g{x, 



C2iv) = ^Tir / dx I def-+-{x,e,v)g{x,e) 



00 Jo 

/ dx / def++-{x,e,v)gix, 



Furthermore, 



9 ■ 2 

+ sm ( 



y/x'^ +v ^ cos^ d + sin^ ( 
x^ + cos^ 9 
+ cos2 9 + w"^ sin^ ( 



is the propagator inverse and 

[ax^ + h cos^ 9 -I- czJ"^ sin^ ( 



r^^(x,0,«) 



(x2 + cos2 9 + "2 sin^ 9)^ 



(A5) 



(A6) 



(A4) 



with a,b,c — ±. This is the full set of equations re- 
quired for the calculation of all running parameters used 
in Sec. [mP] 



Appendix B: CONTRIBUTION OF THE 
EFFECTIVE BOSONIC MODE TO THE HEAT 
CURRENT 

The derivation of the heat current operator in interact- 
ing electronic systems is a long-standing problem stem- 
ming from the fact that a temperature gradient cannot be 
represented as a mechanical perturbation^i^i^^i^. Com- 
monly, a Lagrangian approach is adopted to derive the 
appropriate relations. In the case of an electronic system 
interacting with a bosonic system, the microscopic ex- 
pression for the heat current was derived for the case of 
phonons by Vilenkin et al. —. This procedure, however, 
docs not exactly apply to our problem, due to the fact 
that the effective low-energy action of the bosonic mode 
is created by the electrons themselves. In general, for 
the system described by Eqs. (PT|) . ((^ . and ([^ it is 
practically impossible to disentangle the contributions to 
the energy current due to the electronic and the bosonic 
degrees of freedom. The full action of our problem at 
hand is given by the electronic d-wave superconductor 
coupled to the bosonic nematic mode, both of which can 
carry heat current. In a generic electron-phonon system 
the complete expression for the heat operator was given 
by Vilenkin et al.^^ . Our problem, however, is different 
in the sense, that the effective dynamics of the bosonic 
mode is created by the electrons, see Eq. (|2.1ip . This is 
very similar to problems studied in the context of slave 
particle theories for the t- J model or in composite fermion 
theories of the fractional quantum hall effect, where an 
effective action for a U(l)-gauge field is generated and 
the analysis of transport quantities, especially in the d.c. 
limit, is tedious, see Refs. ,48 ^49..5Q.5L 54 

In order to derive an expression for the effective heat 
current carried by the bosonic mode we expand the sad- 
dle point action to third order following a minimal cou- 
pling scheme in the spirit of the minimal coupling to the 
electromagnetic field. The role of the electromagnetic 
field is in this case assumed by a ther mal gau ge field and 
is explained in great detail in Refs. [30ll32ll33i In the fol- 
lowing we assume a thermal gradient in the direction of 
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the Fermi- velocity at node 1, without loss of generality. 
This implies we have a thermal gauge field, denoted A^, 
without a y-component in the local basis defined at node 
1. Furthermore, we will later make the approximation 



that V ^ 1, which implies that the thermal gradient de- 
couples from the second term in Eq. (|2.ip . We thus find 
the minimally coupled version of the problem to be given 
by 



S = 



I WAfcy<5k,k' dr 



k.k' 



VF{kxSk,k' - iA^{k - k')i d r) 



dr - VpkySuk' 



kit. V ^A(fcx<5k,k' - «Ai(k - k')z d r) 



WA(fc2;<5k,k' - »^i(k - k')i d r) 

dr + VpkySk^ki 



2k'CT > 



r 



Ik', 



(Bl) 



where d r = yd t — d rj, see Ref. |32|. We can proceed 

with the derivation of the effective heat vertex of the 
bosonic mode, which is achieved calculating the diagram 
shown in Fig. [5] a.). The "heat-vertex" with incoming 



frequency f2„ and the external momentum q of the gauge 
field equal to zero, i.e. q = 0, equates in the limit of zero 
temperature to 



r(j.„,k, f]„,q = 0) 



lim 



vNfk. 



1617„ 

Nfkx 

1617„ 



v^kl + v\kl + vl 



\kl + vlkl + {Vn + ^uY 



(B2) 



The derivation of this quantity is very tedious and in- 
volves the use of Feyman parameters, see e.g. Ref. [isl . 
However, it is interesting to note that the above expres- 
sion, before taking the limit 1/v — > in the last line, 
is exact. With the effective "heat-vertex" at hand we 
can proceed to calculate the bosonic contribution to the 
thermal transport. The corresponding diagrammatic ex- 
pressions assume the form shown in Fig. [5]b.). At this 
stage it is very important to note that the above expres- 
sion is even under Vn ^ f2„ — > and k — > — k. 
This is very important, since it implies that the heat- 



current-heat-current correlation function is given by 

r(z/„, k, iin, 0) (r(t/„, k, nn,o) + r(-M„, -k, -ri„, o)) 

G0(k, f„ + r2„)G0(k, Un) 

r(l^„, k, Qn, O fG4,ik, + f^n)G^(k, Vn) . (B3) 

One comment has to be made at this stage, which is, 
that the heat-vertex in our case has been evaluated us- 
ing the free Bogoliubov quasiparticle propagators. This, 
in general, constitutes a problem, once the d.c. limit 
in a response function is taken, since the resulting in- 
tegrals are ill-defined and lead to infinite response coef- 
ficients^^. This is also reflected in the large-N calcula- 
tion in a fermion-gauge field system as they often appear 
in gauge theory descriptions of strongly interacting elec- 
tronic systems of Ref. 150j, where the calculation becomes 
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invalid in the low-frequency limit and the temperature 
has to serve as a cutoff. Starting from a heat-vertex cal- 
culated from free fermionic propagators, a calculation of 
the zero-frequency limit of the bosonic contribution to 
the thermal conductivity using the appropriate propa- 
gator for iJ ^ 1 introduced in Eq. (|2.13p following the 
prescription given in Ref. 'ss'does not yield a finite value. 
This can be traced back to the contribution of the free 
fermions and is an artifact of the neglect of the self- 
energies^. We have not explicitly performed the cal- 
culation taking into account the self-energies. It is, how- 
ever, straightforward using a spectral representation of 
the electronic contributions and introducing a thermal 
broadening^. The important point is that the contribu- 
tion of the bosonic modes to the thermal conductivity 
is finite and naively proportional to , which can easily 
be checked by taking the anisotropic bosonic propagator 
(Eq. (|2.13p ) and scaling all the momenta such as to make 
the resulting integral dimensionless. In this sense the 
bosonic contribution is down by 1 /v with respect to the 
fermionic contribution steming from the nodal fermions 
which move in perpendicular direction to the Fermi sur- 
face, compare Eq. (|3.17p in Sec. IIII Bl However, it is un- 
clear whether this power-counting argument applies once 
the self-energies are taken into account. This question is 
postponed to later works. For a perfectly isotropic sit- 




Figure 5: the effective heat vertex proportional to v according 
to Eq. (|B2[I is shown in a.); b.) shows the bosonic contribution 
to the thermal current. 

nation an analysis of the transport properties in terms 
of the fermionic model seems more appropriate, see Ap- 
pendix [Cl leading to an infinite thermal conductivity in 
the absence of disorder. 



Appendix C: COMMENTS ON FINITE 
CONDUCTIVITIES IN IDEAL SYSTEMS 

Naively, in a clean system neglecting Umklapp scatter- 
ing one would expect infinite response of the system to 
small perturbations, i.e. infinite transport coefficients, if 
interactions are momentum-conserving. 

If we consider a clean Fermi liquid with electron- 
electron interactions we observe the following. If we ap- 
ply an electrical field across the system and forbid Umk- 
lapp scattering, no current is relaxed and the response 



is infinite. On the level of a Boltzmann equation, this 
can be seen very naturally. Assuming electronic quasi- 
particles, we can write the generic linear response form 
of the distribution function as 

/(k)-/0(fc) + /"(fc)(l-/"(fc))Xk (Cl) 

Assuming a generic electron-electron interaction, we can 
write down the following schematic form.— 

X(k) - J dkldqFik, kl,q) [Xk + Ak, - Ak-q - Ak,+q] 

(C2) 

for the scattering integral. It is obvious from the above 
equation that the right hand side vanishes for a particular 
choice of Ak , namely 

Ak = ck , (C3) 

which just restates the conservation of momentum. The 
associated mode is a so-called zero-mode of the scattering 
operator, which renders the matrix inversion singular. 
A zero mode is a mode, which cannot decay, i.e. the 
associated scattering time diverges, which also implies 
the conductivity as defined by the Drudc^^ formula to 
diverge. 

One would expect the same kind of reasoning to ap- 
ply for the thermal conductivity. However, the thermal 
conductivity usually is protected from this divergence 
(Refs. I32II41II4I ) . This can be related to the boundary 
condition of a vanishing electrical current, which forbids 
to excite the corresponding mode. 

It is instructive to compare the above reasoning with a 
relativistically invariant (or equivalently strictly particle- 
hole symmetric) electronic theory, such as the Dirac the- 
ory in the way it applies to intrinsic graphene in the 
clean limit. It was pointed out in Refs. [24I25I that the 
conductivity in such a system can be finite, i.e. a mo- 
mentum conserving interaction can relax a current due 
to the special particle-hole structure of the Dirac Hamil- 
tonian. However, this is not true for the thermal conduc- 
tivity, which is infinite. This can be traced back to the 
conservation of the energy component of the momentum- 
energy tensor and it is a generic property of relativisti- 
cally invariant theories without any sort of translational 
symmetry breaking-!^. Another way to state it is that 
in a relativistically invariant system the current and the 
heat current operator are orthogonal to each other, due 
to particle-hole symmetry, thus no boundary condition 
can cure the divergence. 

This again can very nicely be seen in a Boltzmann 
transport approach. Here again, the collision integral as- 
sumes the generic form shown in Eq. (jC2[l and the exis- 
tence of a zero-mode implies a diverging scattering time. 

In our approach we reformulate a problem of inter- 
acting electrons in terms of electrons interacting with 
bosonic degrees of freedom. However, we will try to re- 
formulate the present problem in terms of fermionic vari- 
ables, yielding a faithful low-energy description. In order 
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to do so we integrate out the effective bosonic degrees 
of freedom. This leads to a long-range electron-electron 
interaction, where the scattering operator for w = 1 looks 
identical to the scattering operator considered in a recent 
publication in the context of intrinsic graphene^^^. In 
this case the thermal conductivity is infinite due to the 
presence of a momentum (energy) mode, which is not 



relaxed due to the presence of electron-electron interac- 
tion. We will show below that this changes in the case 
of anisotropic velocities. In order to do so we consider 
the ansatz introduced in Eq. p. lip , where we chose a 
constant function 'i>p{k,9,v) — ^'^(fc,6',w) — const. We 
can write the linearized version of the scattering integral 
as 



^ « II [v^VT(v^ • k) + v^VT(v^ • k') - v^Vr(v^ • (k + q)) - v^VT(v|, • (k' - q)) 

b 

+ v^Vr(vl-k)+viVr(vi-k')-vlVT(v^.(k + q))-viVT(vi.(k'-q))] . (C4) 



In the case of no anisotropy, this expression vanishes iden- 
tically under local rescaling of the corresponding vari- 
ables. However, if ^ 1 this chances. The contribution 
from the scattering of electrons from one node to the 
same node, i.e. a to a or 6 to 6 can always be annihi- 
lated. However, scattering between adjacent nodes can- 
not, as can readily be seen from analyzing the expression 
corresponding to such an event 

I oc v^Vr(v^ ■ q) (iJ-2 _ 1) 

+ vlVr(vl • q) - 1) . (C5) 

This implies that the energy mode is a zero mode of the 
problem for the isotropic system, i.e. v = 1, whereas 



the mode becomes massive once v'^ ^ 1. This statement 
is equivalent to saying that the interaction with the ne- 
matic mode breaks the relativistic invariance of the un- 
derlying theory, which allows to have non-infinite ther- 
mal conductivity, since the divergence is not protected 
by symmetry any more. We thus expect the scattering of 
electrons from bosonic degrees of freedom (Eq. ()3.28|) 1 to 
correctly describe this scattering for the case of extreme 
anisotropy, i.e. w S> 1, since then the bosonic theory 
seems to be the natural choice. The analysis of a clean 
electronic system and the interaction mediated thermal 
conductivity was carried out in Refs. l40pll and redone 
in a decoupling scheme in bosonic fields by Catelani and 
Aleiner— . 
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